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ABSTRACT 

We consider two types of mechanical coupling between the crust and the core of a precessing neutron star. 
First, we find that a hydromagnetic (MHD) coupling between the crust and the core strongly modifies the star's 
precessional modes when ta <^ (T, x TpY^^; here ta is the Alfven crossing timescale, and and Tp are the 
star's spin and precession periods, respectively. We argue that in a precessing pulsar PSR B 1828-1 1 the restoring 
MHD stress prevents a free wobble of the crust relative to the non-precessing core. Instead, the crust and the 
proton-electron plasma in the core must precess in unison, and their combined ellipticity determines the period 
of precession. Link has recently shown that the neutron superfluid vortices in the core of PSR B 1828-1 1 cannot 
be pinned to the plasma; he has also argued that this lack of pinning is expected if the proton Fermi liquid in the 
core is type-I superconductor In this case, the neutron superfluid is dynamically decoupled from the precessing 
motion. The pulsar's precession decays due to the mutual friction between the neutron superfluid and the plasma 
in the core. The decay is expected to occur over tens to hundreds of precession periods and may be measurable 
over a human lifetime. Such a measurement would provide information about the strong n-p interaction in the 
neutron-star core. 

Second, we consider the effect of gravitomagnetic coupling between the neutron superfluid in the core and the 
rest of the star and show that this coupling changes the rate of precession by about 10%. The general formalism 
developed in this paper may be useful for other applications. 

Subject headings: neutron stars 



L INTRODUCTION 

The most conclusive evidence for a free precession of an 
isolated pulsar comes from Stairs, Lyne, and Shemar (2000, 
SLS); see also Cordes 1993 and Shabanova, Lyne, and Urama 
2001. Their discovery has shown convincingly that some pul- 
sars are precessing, and has opened a new window into the in- 
terior of neutron stars (Link and Epstein 2001, Jones and An- 
dersson 2001, Link and Cutler 2002, Cutler, Ushomirsky and 
Link 2003, Wasserman 2003, Link 2003; see also Link 2002 
for a review). The pulsar PSR B 1828-1 1, which has been mon- 
itored by SLS for about a decade, is spinning with the period 
of 0.4 seconds and precessing with the period of 500 or 1000 
days. The large ratio of the precession to spin periods is dif- 
ficult to reconcile with the current theoretical ideas about the 
neutron star's internal structure. In particular, it has long been 
argued that the neutron superfluid vortices are pinned to the 
crystal lattice of the crust; this has been used to explain pul- 
sar glitches (sudden spin-ups of young isolated pulsars). How- 
ever, as was shown in the pioneering work of Shaham (1977), 
the crustal pinning leads to rather short precession periods. 



t precession " 



■ (^star /^superfluid) ^spin- Here Tspin and Tprecession 



spin and the precession periods respectively, /star and /superfluid 
are the moments of inertia of the star and the pinned superfluid 
respectively. The expected precession period of PSR B 1 828-1 1 
would be of the order of 100 seconds, in sharp contrast with 
what has been observed. Link and Cutler (2001) have proposed 
a way out of this contradiction: they argue that the observed 
precession is so strong that the superfluid vortices are unpinned 
from and move freely through the crustal lattice. Another pos- 
sibility is that the vortices in the superfluid in the crust are never 
pinned to the crustal lattice; this has been argued on theoretical 
grounds by Jones (1998). We feel that both of these ideas, while 



potentially viable, require more detailed calculations. 

In addition, if the proton-electron plasma in the core par- 
ticipate in the precessing motion and if, as is commonly 
believed, the protons condense into type-II superconductor 
(Baym, Pethick, and Pines 1969), then the expected strong in- 
teraction between the superconductor's fluxtubes and the vor- 
tices of the neutron superfluid does not allow slow precession 
with small damping (Link 2003). The conflict with observa- 
tions is avoided if either : 1 . the proton-electron plasma does 
not participate in the precessing motion, and the crust alone 
precesses ("Chandler wobble"); we will show that this possibil- 
ity is excluded due to the MHD crust-core coupling, or 2. the 
protons in the core do not form type-II superconductor, as is 
commonly believed, but instead they form a type-I supercon- 
ductor (Link 2003). This is not far-fetched since proton pairing 
calculations in the core are uncertain; moreover, recent work 
by Buckley, Petlitski and Zhitnitsky (2003) argues that the in- 
teraction between the proton and neutron condensates may turn 
the neutron-star interior into type-I superconductor even if the 
proton pairing calculations favor the type-II phase. 

While it would be exciting to gain some observational han- 
dle on the exotic quantum fluids in the neutron-star interior, the 
main focus of this paper is elsewhere. Here, we concentrate in- 
stead on the MHD and gravitomagnetic coupling between the 
crust and the core; these effects make an impact on the preces- 
sion dynamics yet their nature is well-understood theoretically 
and (in case of MHD) is well-tested in laboratory experiments. 
We will, however, also discuss the decay of pulsar precession 
due to the mutual friction between the neutron superfluid and 
the proton-electron plasma in the core. 

This paper is structured as follows. In section II we present 
a toy model for MHD coupling between the crust and the core 
and solve the Euler's precession equations within this model. 
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We show that the precession period is strongly affected once 
the timescale of magnetic coupling is comparable to the geo- 
metric mean between the spin and precession periods. In this 
section we also consider the damping of precession via elec- 
tron scattering off the magnetized neutron superfluid vortices 
in the core (this precess is called mutual friction). The damping 
timescale is generally tens or hundreds of precession periods, 
and its exact value is sensitive to the effective proton mass in 
the core. Thus by monitoring the precession decay over a few 
decades one may be able to constrain the strong n-p interactions 
in the neutron-star core. 

In Section III we move away from toy model and consider the 
nature of slow MHD waves in a rotating gravitationally strati- 
fied neutron-star core. Our calculations for this more realistic 
model generally confirm our toy-model results. 

Finally, in section IV we take into account relativistic frame- 
dragging around spinning neutron star, and find two modes of 
relativistic precession. The first mode is the Lense-Thirring 
(LT) precession of the crust in the gravito-magnetic field of 
the core, first considered by Blandford and Coppi (Blandford 
1995). The Lense-Thirring precession is relatively fast; its pe- 
riod is only an order of magnitude greater than the spin period 
of the pulsar. The LT precession is most easily excited when 
the crust suddenly changes its angular momentum (e.g. due to 
a collision with an asteroid). It is damped on the timescale of 
viscous coupling between the crust and the core, about 100 sec- 
onds. 

The second mode is the Eulerian precession which is mod- 
ified by the inertial frame dragging. This mode is most easily 
excited by a change in the crustal tensor of inertia, e.g. by a 
sudden deformation of the crust due to magnetic forces. We 
find that the frame dragging modifies the precession frequency 
by about 10%. 

2. IDEALIZED MODEL FOR THE CRUST-CORE MHD COUPLING. 



Ohmic dissipation inside the neutron star is very slow com- 
pared to the precession period, and therefore ideal MHD pro- 
vides an excellent description of the neutron-star interior. The 
magnetic field threads both the crust and the core; in an ideal 
MHD the relative displacement of the crust and the core cre- 
ates magnetic stresses which oppose this displacement. Thus, 
the magnetic field lines act as elastic strings (e.g., Blandford 
and Thorne 2004). Motivated by this we follow the spirit of 
Bondi and Gold's (1955, BG) analysis of the Chandler Wobble, 
and consider the crust and the core as solid bodies coupled by a 
torque which opposes their relative displacement: 



(1) 



Here S4> is the small angular displacement between the crust 
and the core, and /i is a constant representing the strength of 
MHD coupling. While this model is simplistic, it is (a) fully 
solvable and (b) correct in predicting the main features of the 
precession with MHD crust-core coupling. We consider a more 
realistic model in the next section. 

For mathematical simplicity, we assume the crust is axisym- 
metric and we work in the coordinate system (e^i ,62,^3) rigidly 
attached to the crust so that 63 is directed along the symmetry 
axis. We also assume the core to be spherically symmetric. We 
denote by (A,A,C) and {D,D,D) the crust's and the core's three 
principal moments of inertia, respectively. The dynamics of the 
system is described by the coupled Euler's equations [cf. Eq. (1) 



ofBG]: 

Aa}i + (C—A)lu2W3 = —^5(j)\ = —D{fli +Ci;2^^3~ ^3^^2)1 
Auj2 — (C — A)luiU!3 = —fi64>2 = —D(fl2 — LUiflT, + uj3Qi), (2) 

Cl03 = -fl54>3 = -DCI3. 

Here uj = (uji,uj2,(^3) and (l = (rii,rj2,f^3) are the angular ve- 
locity vectors of the crust and the core, respectively, and the 
sign of d<j) is chosen so that 



dScf) 
dt 



= uj — fl. 



(3) 



The observed wobble angle of PSR B1828-11 is only ^ 3 de- 
grees (Link and Epstein 2001); motivated by this we restrict 
our analysis to small-amplitude precession. More precisely, we 
consider a small periodic perturbation of an equilibrium state 
in which both the crust and the core are rotating around the 
crust's symmetry axis with the angular velocity n. The dynam- 
ical quantities are then expressed as follows: 



and analogously. 











L02 




L03 


= n + uj3e 






^2 


= n2e'-', 



(4) 



(5) 



n3 = n + n3e'' 



Here it is assumed that the complex amplitudes w,- and 17,- are 
small compared to n. It is also convenient to define, in the usual 
way, the crust's ellipticity: 



C-A 



(6) 



In the dynamical Equations (|3 we can neglect the terms which 
are of second order with respect to a), and f2,, and use Eq. 
to eliminate d(j>. The linearized equations of motion are given 
below: 

/i ~ 

iaACj\ + eAnuj2 = — (Q,i-uji) = —iDaVL \ —Dn{u}2— ^i), 
ia 

iaAu!2-eAnCji = —(Cl2-Lij2) =-iD(TCl2 + Dniu!\-Cl\), (7) 
ia 



iaCCj3 = —(CI3-UJ3) ■■ 
ia 



-iDafl 



The third equation above is decoupled from the first two; it de- 
scribes the small rotations of the crust and the core around the 
symmetry axis of the crust. This equation alone gives two fre- 
quency eigenvalues: 



^3 = 0, 



(74 = 



H(C + D) 
CD 



1/2 



(8) 



The trivial eigenvalue 0-3 corresponds to the crust and the core 
rotating in unison without any relative displacement, whereas 
the eigenvalue CT4 corresponds to the crust-core oscillations 
around the rotation axis; these oscillations are not affected by 
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the ellipticity of the crust and the rate of stellar rotation, and 
do not represent precession. The information about preces- 
sion is contained in the first two equations of 0. Follow- 
ing BG, we can simplify the algebra by considering the sum 
(first equation)+;x (second equation), and by introducing the 
new variables, = uji +iCj2 and ri"*" = ^li +iVl2. In the end, 
we get the following eigenvalue equation: 



(a-m){a +na) = fi 

This equation has three solutions: 

(To ~-n-aljn, 

0'l,2 



A + D n 



(9) 



/2, 



where 



/^(A + D) 



AD 



(10) 
(11) 

(12) 



is the frequency of the mode in which the crust and the core 
of a non-rotating star oscillate differentially, with the restoring 
force of purely MHD origin, and 



(13) 



Since in our case A <C Z), one can show that ap ^ 2cro for all 
values of e and /i. We therefore have 




CTi ~ (Tp = en + al/n, 



A+D 



en(l + en^/f7^,) 



(14) 
(15) 



The frequency cti characterizes the differential precession be- 
tween the crust and the core; in the limit of zero magnetic cou- 
pling (i.e., a,„ = 0) its value tJi = ne is the frequency of a free 
precession of the crust. By contrast, the frequency (72 corre- 
sponds to the mode in which the crust and the core are trying 
to precess in unison. In the limit of infinite magnetic coupling 
(i.e., (T,„ = oo) its value of (72 = enA/(A + D) is the frequency of 
precession of the neutron star as a whole: the crust is the source 
of ellipticity but the core is rigidly attached to the crust and they 
precess together. 

Let's apply these results to PSR B 1828-11. The infeiTed 
dipole magnetic field of this pulsar is Z? ~ 5 x lO'^G (see SLS), 
and the Alfven speed in the core is 



B 



2 X lO'Vcm^ 



1/2 



5 X IQi^G, 



cm/s (16) 



when the core is not superconducting, and 



1.4 X 10" 



B 



5 X 1012G 



1/2 



2 X lO'Vcm' 
P 



1/2 



cm/s 



(17) 



when the core is superconducting. Here p is the density of the 
core material which is participating in the Alfven-wave mo- 
tion. We estimate the characteristic cr,,, ^ ttVo/R to be 0.3s"' 
for a non-superconducting core, and 4.2s"' for a superconduct- 
ing core. In deriving these numbers we have assumed that all of 
the core is participating in the Alfven-wave motion; we remark 
that if the neutrons form a superfluid, then only the charged 
proton-electron plasma is magnetically coupled to the crust, and 
the estimates for (j„, should increase by a factor of ^ 4. The 
spin period of PSR B 1828-1 1 is T^pin ~ 0.4s, and from Eq. ^ 
we see that the period of the crust-core differential precession 
(Chandler Wobble) is 



Ti =27r/(7i - lO^^s 



for a non-superconducting core, and 



Ti -5s 



(18) 



(19) 



for a superconducting core. In the above estimates, we have 
assumed zero ellipticity for the star and that all of the core is 
magnetically coupled to the crust; thus our estimates are upper 
limits on T\. Since the precession period of PSR B1828-11 is 
— 4 X lO^s, we can say with certainty that the observed preces- 
sion is not the "Chandler wobble" of the crust relative to the 
core. Rather, in agreement with the argument sketched by Link 
(2003), the crust and the magnetically-coupled part of the core 
precess in unison' and their precession period is found from 
Eq. IB): 



Tspin ^+^_o„in7( ^spin 



10"^ \ f\+D/A\ 

(20) 

In deriving Equations M5\ and (I20> . we have assumed that there 
are no extra torques acting on the charged plasma of the core. 
This assumption breaks down if the core is a type-II supercon- 
ductor and the neutron superfluid vortices interact strongly with 
the magnetic fluxtubes. Link (2003) has shown that a strong 
vortex-fluxtube interaction is inconsistent with the observed 
precession on PSR B 1828-1 1, but has pointed out that the dif- 
ficulty is alleviated if the core superconductivity is of type-I 
rather than type-II. In this case, the magnetic field is contained 
not in quantized fluxtubes but in larger domains (although these 
domains probably still thread densely the neutron-star interior). 
Then, the relative motion of the plasma and the neutron super- 
fluid is damped by the scattering of the electrons on the magne- 
tized supefluid vortices [Alpar, Langer, and Sauls 1984, Alpar 
and Sauls 1988]. This damping is known as "mutual friction", 
and its characteristic timescale is 

fmf = lOT^pin(mp/Sm*ff{m*,m*,A„,pc,p), (21) 

cf. Eq. (32) of Alpar, Langer, and Sauls. Here mp/m„ and ni*/m* 
are the bare and the effective proton/neutron masses, A„ is the 
neutron condensate gap, pc is the density of the proton-electron 
plasma in the core, and / is a function which depends weakly 
on its variables. Alpar and Sauls (1988) give detailed discus- 
sion of fmf, and we refer to them for the details. 



'Link (2003) has argued that the crust and the charged part of the core precess in unison when the precession frequency is o^prec < cr,,,. However, this estimate does 
not talce into account the rotation of the star; we see from our Eq. 1141 that the correct criterion is ojprec < fm/"' ^ rnore stringent condition. 
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The precession damping timescale Tp, 
following relation^ (see, e.g., BG): 



is determined by the 



Tpr = 72fmf/7]ipm ^ I5(mp / Sm*)^yeeLTS. 



(22) 



Thus, the mutual-friction damping of precession may be ob- 
servable for PSR B 1828- 11 over the timescale of human life, 
and its measurement will yield the information about Sm*/mp. 

3. ALFVEN WAVES IN A ROTATING NEUTRON STAR. 

It is interesting to note that in the absence of the crust ellip- 
ticity, the frequency of the neutron-star Chandler Wobble scales 
as 1/n, see Eq. (I14> . As is seen from this equation the fast ro- 
tation reduces the effectiveness of MHD crust-core coupling. 
For a neutron star with a fluid core the coupling is mediated 
by the Alfven waves which are excited by the precessing crust 
and propagate into the core. The characteristic timescale for 
the coupling is ^ R/vap, where Vap is the speed of these Alfven 
waves. Thus, we expect that the Alfven waves are slowed down 
as the star spins faster; our detailed analysis below confirms this 
expectation. 

MHD in rotating fluids has been the subject of an extensive 
research in geophysical fluid dynamics, with applications to the 
Earth's fluid core [see Hide, Boggs, and Dickey (2000) and ref- 
erences therein.] There is, however, a significant difference be- 
tween the Earth and neutron star cores. The Earth interior is 
approximately isenthropic, and the Taylor-Proudman theorem 
is applicable; thus the velocity field is almost constant along the 
lines parallel to the rotation axis. By contrast, the neutron-star 
interior is stable stratified due to the core's radial composition 
gradient (Reisenegger and Goldreich, 1992). The fluid motion 
is restricted to equipotential shells, which we assumed to be 
spherical (this is a good approximation for the slowly-spinning 
PSR B 1828- 11). The motion is strongly subsonic, therefore the 
velocity field is divergence-free. 

Under these conditions, the general small fluid displace- 
ments can be represented by the radius-dependent stream func- 
tion ip(r, 9, (p), so that in spherical coordinates the displacement 
components are 



C = 0, 



(23) 



1 



The radial component of the vorticity of the fluid is 

7? = (i/r)2(a/a)v2v, 

where is the Laplacian operator on the unit sphere: 



(24) 



radial component of the absolute vorticity, cf. Eq. (5) of Levin 
and Ushomirsky (2001): 



dt 



(77+2^1 cos 0) = r-W xds. 



(26) 



Here d/dt = d/dt + v-W is the Lagrangian time derivative, r is 
the unit radial vector and gb is the acceleration due to the restor- 
ing magnetic stress. Following Kinney and Mendell (2002), we 
restrict ourselves to the special case of the spherically symmet- 
ric radial magnetic field, B = B(r)r. The results obtained below 
should be qualitatively correct for a more general field configu- 
ration; however, we have chosen a particularly simple geometry 
in which the mathematical evaluation of the right-hand side in 
Eq. ( I26> is greatly simplified. The relevant components of dg 
are given by 



Ob ■ ee.i 



1 d 
Anpr dr 



d_fCi 
dr 



(27) 



We can now write down the linearized equation of motion for 
the stream function: 



52 
dt 



B'r' 



dr 



(28) 



We look for the solution of Eq. (|28 

7/.(r,0,0) = I], 

Since Yi,„ is an eigenfunction of both d/d^) and V^, Eq. ( I28> 
separates into individual ordinary differential equations for V'/mj: 



in the following form: 

^Ur)YU0,c[>)e'''>"'. (29) 



2mVLa, 



Im 



>lm~ 



1 d 
Vim{r)+- — 

47r/3 or 



bV 



d_ ( i^imir) 
dr 



= 0. 
(30) 

We now consider the short-wavelength (WKB) approximation 
for the above equation, and hence derive the following disper- 
sion relation: 



Im 



(31) 



1(1+1) 

where k is the radial wavevector The purely toroidal Alfven 
waves correspond to the case when m = in the above equa- 
tion. These waves are not affected by the stellar rotation and 
have the dispersion relation identical to that of Alfven waves 
in a non-rotating star. However, in PSR B 1828- 11 the Alfven 
waves are excited by the slowly precessing crust, and therefore 
one should consider the waves with / = 1 and m = -l. In this 
case the second term on the right-hand side of Eq. ( 13 1> is the 
dominant one since a fl, and the wave is strongly slowed 
down by the stellar rotation, just as we expected. The radial 
wavelength of the excited Alfven mode is 



Aa = 27r/^~27rVa/A 



(32) 



V = 



1 



smt 



d_ 

d9 



sin9- 



dO 



smt 



(25) 



which equals ^ 6 x 10 cm and ^ 8 x 10 cm for normal and 
for superconducting neutron star interior, respectively. In both 
cases, it is more than two orders of magnitude greater than the 
radius of neutron star, ^ 10''cm. Therefore, the part of the 
neutron-star interior which is magnetically coupled to the solid 
crust will precess in unison with the crust. This conclusion is 
robust and is in agreement with our toy-model results from the 
previous section. 

^Alpar and Sauls (1988) have erroneously overestimated the precession damping timescale by a factor p/ pc- They have associated the viscous damping timescale 
with txaipl pc, since this is the timescale it takes for the neutron superfluid to come to co-rotation with the charged plasma. However, the neutron superfluid carries 
most of the star's moment of inertia, and if the superfluid spins at a different rate than the rest of the star, it is the crust+plasma which are coming to co-rotation with 
it. Therefore one should use for the viscous damping timescale when one evaluates the precession damping timescale. 



Since the fluid inside the neutron star is strongly stratified by 
gravity (i.e., the Brunt- Vaissalla frequency ^ S^jC,,,), the 
fluid motion on different shells is coupled only via magnetic 
stresses. One can write down the dynamical equation for the 
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4. GRAVITOMAGNETIC COUPLING BETWEEN THE CRUST AND 
THE CORE. 

The gravitational redshift at a neutron-star surface is ~ 0.3, 
and therefore relativistic effects, including the dragging of the 
inertial frames, are strong in and around neutron stars. In this 
section we analyze how the frame-dragging affects the relative 
precession of the crust and the core. Our post-Newtonian cal- 
culations rely on the usage of the gravitomagnetic field, H, see 
Thorne, Price, and MacDonald (1986) for the details of this for- 
malism. 

4. 1 . The gravitomagnetic coupling torque 

Consider the gravitomagnetic force acting on the small re- 
gion of the crust of mass dm^ . In the post-Newtonian approxi- 
mation, it is given by 

dFcM = dmvx H = dm{uj xr)xH, (33) 

where v*, H, uj, rare the velocity of the small region, the gravit- 
omagnetic field, the instantaneous angular velocity of the crust, 
and the radius-vector of the region, respectively. The torque 
acting on this region of the crust is given by 

dT = Fx dFcM = dm(H ■ r)uj x r, (34) 

where we have used the vector identity {AxB)xC = {C ■ A)B - 
(C ■ B)A. The field H is that of a dipole, and 

H ■ r = -(4/r^)/- r = -(4D/r^)Q ■ r, (35) 

where J, d, and D are the angular momentum, the angular ve- 
locity, and the moment of inertia of the spherical core [Here 
we ignore interaction of the crust with its own gravitomagnetic 
field. It can be shown (Thorne and Gursel, 1983) that this self- 
interaction can be absorbed into the free precession.] 

Now, substituting Eq. i35i into Eq. ( I34> . and integrating over 
the crust, we arrive to the following form of the gravitomagnetic 
torque: 

fgn, = a;x/gn,r2, (36) 

where is the linear operator (represented, generally, by a 
3x3 matrix) defined as follows: 

I„^n = - j £/Vp(r)(4D/r^)(l1 • (37) 

We use Dirac's the bra and ket notation and express this opera- 
tor as 

h^ = - j d\p{r)(ADlr')\f><f\. (38) 

From the above expression we see that /g^ is a hermitian op- 
erator: since the integrand cx |r >< r| in Eq. ( I38l l is hermitian, 
the integral must also be hermitian. This means that the matrix 
representing /gn, is symmetric. 

If the crust is spherically symmetric, then /gn, = pi, where p 
is a real number and / is a unit matrix. In this case, the torque 
acting on the crust is 

f = pu;xn, (39) 

which is the familiar form of the Lense-Thirring torque acting 
on the gyroscope. In the situation considered here the crust is 
slightly deformed, so that 

Ig^ = pI+epK, (40) 

where /iT is a 3 x 3 matrix with entries of order 1 . 

'in this section we shall collectively refer to the solid crust and the core plas 



4.2. The dynamics of relativistic precession 

Again, we use the BG approach to the precession of an in- 
teracting crust and core. The equations of motion which in- 
clude the gravitomagnetic torque components Ti,T2,Tj, are [cf. 
Eq. (l)ofBG]: 

Alui + (C — A)u!2LiJi = X(Qi —LUi) + Ti = —D[Qi + W2^^3 — ^3^2], 

Auj2 - (C-A)uJiUJi = A(f72 - UJ2) + T2= -D[n2 + ^3^1 - wi f^#) 

Cw's = X(flj, - uj^) + T^ =-D[flj + u!ifl2-uJ2^i]- 

Here we have added the terms A(f2 - uj) which represent the 
viscous torque between the crust and the core (e.g., due to mu- 
tual friction between the neutron superfluid and the core plasma 
coupled to the solid crust). As in the previous sections, we are 
interested in the small amplitude precession, when the motion 
differs only slightly from the rigid motion rotation about z-axis, 
so that oji, UJ2, UJ3 = uj3-n, rti, ^2, and (I3 = r23-« are all much 
less than n. Then Eq. J41> becomes 

AuJi + (C-A)nuj2 = X(ni-LJi) + Ti=-D[ni+ n{uj2 - ^2)] , 
Aw2-(C-A)«wi = \{n2-uj2) + T2 = -D[VL2 + n{n^-uJiM2) 

CuJi = x(n3- 0)3) + 7^3 = -on^ . 

The general strategy now is to identify the leading terms in 
T[, T2, and Tt,, using Eq. ( I36> . and then solve Eq. J42> . Since 
D A, we consider a simplified case when the spherical core 
has an infinite inertia: D 00, so that the core's spin does not 
change in the inertial frame of reference. Therefore we have, 
from Eq. ( I42> . 

ni + n(uj2 - f^2) = ^^2 - n(wi - f^i ) = fis = 0. (43) 

We look for a mode with the growth rate 7, so that Cti = 
7i7i, etc. By considering the sum [the first component of 
Eq. (l43H -i-ixrthe second component of Eq. ( I43H . we get 

f2+ = -^uj^^ (44) 
in + 7 

where fi""" = ili + 1^2, w"*" = wi + iuj2- 

Let us restrict ourselves to the case of the axially symmetric 
crust. In this case, both tensor of inertia and /g^, diagonalize 
in the same basis because of the axial symmetry. We can then 
write 

/gn,f2 = pniei+pn2e2 + pil + ek)n-ie-i, (45) 

where el 2.3 are the unit vectors along x, y, and z axes respec- 
tively, with the z axis chosen to be the axis of symmetry [as in 
Eq. (I42H . and ^ is a number of order 1. Then the gravitomag- 
netic torque in Eq. ( I36> . to leading order, is 

f = np{[L02{l + ek)-n2] ei + [fii - wi(l + ek)]e2} . (46) 

Now, let us substitute this expression into Eq. ( I42> . add [first 
row]+/ X [second row], and ignore the right-hand side with D in 

coupled to it as the "crust". 
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it [we have already taken care of it by setting D oo]. We 
get, after dividing by A, and substituting 7 instead of the time 
derivative: 

7^;+ + ine[l - {p/A)k'\uj^ = (1 /A)[A + in/7](fi+- uj+). (47) 

Notice that in the above equations the contribution due to the 
gravitomagnetic terms can be represented by an effective mod- 
ification of the elHpticity e e[l -{p/A)k] and of the viscous 
coupUng coefficient A X + inp. Therefore one can consider 
the precession solution without gravitomagnetic terms and then 
substitute the ellipticity and coupling in this solution by their 
modified values. The resulting expressions then represent the 
precession solution which includes the gravitomagnetic cou- 
pling. 

The rernaining calculation is straightforward. Let e = e[l - 
(p/A)k], A = (l/A)[X+inp]. By substituting Eq. (|44j into 
Eq. i4-H . we get the following equation for the growth rate 7: 

-/^+[in(l-e) + X]-/ + n^e = 0, (48) 
which has two solutions: 

7 = ( 1 /2) |-[in( 1 - e) + A] ± \/[i«( 1 - e) + A] 2 - 4n2e I . (49) 

We now use the fact that e <C 1 in Eq. i49l . and we get to the 
leading order in e for the "H-" solution 

H6 - 

T = — 7TTT - (50) 
1 -i(A/n) 

This solution corresponds to the Eulerian precession of the 
crust, with the frequency 

LU,, = ne-m(p/A) = ne[l- {p/A)k] [l-(p/A)], (51) 

and the damping rate 

l/Tp, = (X/A)e[l-(p/A)]. (52) 

The contribution of the frame dragging comes in through terms 
which contain p/A. The frame dragging reduces both the pre- 
cession frequency and the damping rate by relative order of 
p/A. Now, from Eq. ( I39l l, we can work out that p/A = ojur/n, 
whereo^LT is the frequency of the conventional, gyroscopic 
Lense-Thirring precession. The calculations of Blandford and 
Coppi (Blandford 1995) show that uj^j/n ^ 1/7. Therefore, 
we expect the dragging of inertial frames reduces the frequency 
and the damping rate of precession by ^ 10%. 



What about the second, "-" solution of Eq. J48> ? We have, 
to the leading order in e 

ry = -in-X = -in{l+p/A)-X/A. (53) 

This mode corresponds to the situation when the spin of the 
crust is misaligned with that of the core. In the inertial frame of 
reference (as opposed to the frame attached to the body), one 
must take the -in out of 7. The piece that is left is then 

7inertiai = -inp/A-X/A. (54) 

This corresponds to the Lense-Thirring precession considered 
by Blandford and Coppi (Blandford 1995), which is damped on 
a short timescale A/A, i.e. the viscous time on which co-rotation 
of the crust and the core is enforced. 

5. DISCUSSION 

In this paper we have analyzed the effect of the crust-core 
coupling on the Chandler Wobble of the neutron-star crust. 
We have found (section IV) that the gravitomagnetic crust-core 
coupling does not affect strongly the Chandler Wobble, but in- 
stead modifies its frequency by about 10%. By contrast, we 
have found that the MHD interaction between the crust and the 
core of a rotating neutron star dramatically changes the dynam- 
ics of the Wobble, for typical values of the pulsar spin and mag- 
netic field. In particular, we have shown that the observed pre- 
cession in PSR B1828-11 can not be the Chandler Wobble of 
its crust; instead, the crust and the plasma in the core must pre- 
cess in unison. The precession is damped by the mutual friction 
in the core. This damping has a timescale of tens or hundreds 
of precession periods, and may be observed over the span of 
human life. The measurement of the damping timescale would 
constrain the value of dm*/mp and thus provide information 
about the strong p-n interactions in the neutron-star core. 

While the immediate astrophysical impact of our paper is 
modest, it presents some novel analytical techniques. In sec- 
tion III, we have developed the theory of slow Alfven waves 
in a gravitationally stratified uniformly rotating fluid (as is ap- 
plicable for a neutron star). In section IV, we have analyzed 
the precession dynamics of a biaxial rigid body in the presence 
of strong gravitomagnetic field. As far as we are aware, both 
of these technical developments are new, and we envisage their 
further applications to the dynamics of neutron stars. 
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